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Abstract 

In this paper we investigate the effect of introducing transverse mo- 
mentum cutoffs on the BFKL equation. We present solutions in moment 
space for various models of the BFKL kernel for different combinations of 
these cutoffs. We improve on previous calculations by using the full BFKL 
kernel (rather than simplified analytic approximations). The significance 
of the next-to- leading or "higher twist" terms in the kernel are assessed. 
We find that, while these terms are negligible in the absence of cutoffs, 
introducing an infra-red cutoff markedly enhances their significance. 



1 Introduction 

The BFKL equation ]]J can be formulated as an integral equation which deter- 
mines the evolution (in x) of the unintegrated gluon distribution function f(x, k 2 ). 
It is expected to be relevant for sufficiently small x, (i.e. a s \nl/x ~ 1) provid- 
ing the, as yet unquantified, sub-leading corrections are not too large. Cross- 
sections which are sensitive to the small- x dynamics are constructed by convolut- 
ing f(x,k 2 ) with the appropriate coefficient functions (impact factors). In this 
paper we consider the integrated structure function 

rQ 2 dk 2 

j o ^f( x , k 2 ) = G(x,Q 2 ), (1) 

which is equivalent to the usual DGLAP [0 gluon density, xg(x, Q 2 ), in the double 
leading logarithm approximation. However, this is not a physical observable and 
computation, for example, of the deep inelastic structure functions F 2 (x, Q 2 ) and 
Fl(x,Q 2 ) requires less trivial integrals over k 2 . 

Here we do not consider the phenomenologically more interesting cross-sections. 
Rather, we wish to elucidate more on the technical details which are involved in 
computing G(x, Q 2 ) using a modified BFKL kernel. Our modifications corre- 
spond to restricting the virtualities of the t-channel gluons in the real emission 
part of the BFKL kernel, in line with possible constraints imposed by energy 
conservation and the desire to keep away from the region of low virtualities, i.e. 
we impose fixed infra-red and ultra-violet cut-offs on the real emission part of the 
kernel. We introduce these phase-space restrictions on the full BFKL kernel. 

First recall the (angular averaged) BFKL equation for the evolution of the 
(unintegrated) gluon density. It may may be written 



„, , 9 . „ rm , , 3a s k 2 f x o dx' f°° dk' 2 

f(x,k 2 ) = fW( x ,k 2 ) + ^— — — x 

71 Jx X JO k 

Uf(x',k' 2 )-f(x',k 2 )) , f(x',k 2 ) } 



\ \k' 2 -k 2 \ ' Vfc 4 + 4fc' 4 J ' (2) 

Although this equation has been known for some time there are various theoretical 
difficulties in determining its phenomenological implications. The equation is 
derived in a limit where a s is kept fixed, and only includes leading terms (i.e. 
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terms with equal powers of a s and ln(l/x)). Furthermore, the equation involves 
an integration over all transverse momenta including the lowest values where 
perturbation theory is expected to break down (because large radiative corrections 
lead to a s > 1 ). This region may be regulated by introducing an infra-red cutoff 
[1j H> HI- 111 addition, energy conservation demands that an ultra-violet cutoff is 
present ||. The purpose of this paper is to determine the effect that these cutoffs 
have on the mathematical structure of the equation and its solution. 

Clearly in order to produce a realistic model of a small x physical cross section 
some attempt must be made to parameterize the region of phase space which we 
have removed (below the infra-red cutoff) where confinement effects are expected 
to become important. Some phenomenological attempts have been made in this 
direction @, ||. In this paper we restrict ourselves to a consideration of what we 
feel may be reliably calculated within the framework of perturbative QCD and 
make no attempt to discuss what happens outside of this region. 

The equation is best solved by taking Mellin moments with respect to x, k 2 
defined by 



/(*,* 2 ) = / ^r^fix^) (3) 

J c—ioo Z7Y1 

f(x,u) = / dk 2 (k 2 )-i-"f(x,k 2 ) (4) 
Jo 

rc+ioo A AT 

f(x,u) = / x n ^(N 1 uj) (5) 

J c—ioo Z7TZ 

F(N,u) = I dxx- x - N f(x,u). (6) 
Jo 

The contour in eq.(^) lies parallel to the imaginary axis and to the left of 
all singularities in the N— plane. In eq.(|3]) the contour also lies parallel to the 
imaginary axis and to the left of all singularities associated with the small k 2 
behaviour of f(x, k 2 ) but to the right of singularities associated with its large k 2 
behaviour. This choice ensures that the respective inverses of these transforms 
exist and are consistent. The addition of 1/2 to the power of k 2 included in 
eq.(|3|) is merely for convenience. With this definition the tu-plane contour lies 
along dte(uj) = 0. 
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The following solution to the BFKL equation (in 'Double Mellin Transform' 
space) is then obtained: 



where K{u) is the (Mellin transform of the) BFKL kernel 

K(u) = a s [-2 7 - ^(~ +u)- V(i - w)] = a sX {co), (8) 

ip and 7 are the logarithmic derivative of the gamma function and the Euler 
constant respectively and a s = 3a s /7r. The full kernel has an infinite set of simple 
poles at u — { ±1/2, ±3/2, ■ ■ • } as a result of the poles in the ip functions. Fig.fll]) 
shows the real part of the kernel in the region of interest closest to the contour. 

In a previous paper O we presented a method for finding J-(N, u) in the 
presense of infra-red, Qq, and ultra violet, Q\, transverse momentum cutoffs on 
the real terms (i.e. the f(x, k 2 ') terms) in eq.@) for a general input distribution 
FW(N,u). The results may be summarised as follows: 

The function S(N,u) depends in general on Ql and Q\, and upon the pa- 
rameterisation chosen for the input distribution. Explicitly we considered three 
cases: the non-cutoff case, the single (infra-red) cutoff case and the double (infra- 
red and ultra-violet) cutoff case (hereafter denoted with a subscript n, s and d, 
respectively). S s (N,u) is uniquely determined by the fact that it must cancel all 
poles from T S {N 1 u) to the right of the contour in the u;-plane. Similarly Sd{N, u) 
must cancel all poles to the right and to the left of the contour. The structure 
of the solutions in (x, fc 2 )-space depend on the pole structure in Mellin space, the 
positions of which are determined by the zeros of the denominator of eq.(|9]): 

N + K(uj) = 0. (10) 

In this paper we present solutions in (N, fc 2 )-space for various models of the 
kernel and cutoff combinations. Initially we consider a simplified model for the 
kernel, introduced by Collins and Landshoff (§], in which only the two poles 



Im (w) 




Re (K (w) ) 



Re (w) 



Figure 1: Real part of the BFKL kernel for a s = 0.2. Note the saddle point at 
to = (0, 0) and the poles ato; = ±|,±|,---. The poles in T{N, u) in the cj-plane 
are given by the solution to iV + K (u) = 0. 
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nearest to the contour in eq.(^) are kept. These are the poles which lead to the 
dominant behaviour of G(x, Q 2 ) in the limit Q 2 — > oo. We present solutions for 
a particular input distribution, chosen such that its parameters are distinct from 
those subsequently used for the cutoffs on the BFKL evolution. We repeat the 
analysis implementing the exact kernel of eq.®, again keeping only the leading 
tu-plane poles. The integration over the k 2 is performed in iV-moment space 
and the results are presented in terms of the moment of the integrated gluon 
distribution G(N, Q 2 ) || || for each case. We present the effective small-x slope 
in each scenario. 

The inclusion of additional ( "higher twist" ) cu-plane poles in the kernel is then 
implemented (i.e. the two nearest poles on each side of the contour in eq.(|3])). 
Their numerical significance (as a function of x and Q 2 ) is assessed. Again we 
use the full BFKL kernel. 

For the case with no cutoffs present we find that the small-x behaviour, even 
for fairly low values of Q 2 , is dominated by the leading pole in the u-p\ane (the 
next-to-leading poles contribute less than 1% for x < 10~ 3 ). That we can safely 
neglect all poles other than the two which lie nearest the cu-plane contour is 
essential for the validity of the DGLAP formulation of the small- x structure 
functions (e.g. see |, 1, H|, H H) . 



In the presence of infra-red cutoff the addition of the subleading pole (there 
is only one pole to add for an infra-red cutoff) has a more significant effect. We 
present results on the size of this effect as a function of x and Q 2 and comment 
on its possible significance. We find that the next-to-leading terms may be as 
large as ~ 20% for moderate values of Q 2 . 

2 The simple model revisited 

Let us briefly recap the examples given in []7|. The following two pole model for 
the kernel was used 

41n2a s K 
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for which eq.flTUD has only two solutions which we denote inj^i, where 



1 1 ^ K ° 
UN1 = 2 V 1 + ^ 

_ - 2 ~\ 



1 + iV- 1 ^^) UJ 2 - L) 2 



(12) 
(13) 



A I 



In the infra-red cutoff case S S (N, u) removes the pole at u — tom, i-e. it 
satisfies 



J*»(N,u m )+S,(N,u m ) = 0. 
So, the solution is given by 

1 



(14) 



x 



1 + N-iKiu) 

+ (15) 



w- 5 



for k 2 > Qq and zero otherwise. 

In the double cutoff case Sd(N,uj) removes both poles in eq. (|T3|) , i.e. it also 
satisfies 



^°)(JV, -u N1 )+S(N, -u m ) = 0, 
and the solution is given by 



(16) 



1 + N^Kiu) 



x 



^\N,uo) + {A(oo n1 )S(N,lu)} 



for Qq < k 2 < Q\ and zero otherwise, where 



A(u n1 )S(N,uj) 



Q 



-2u) 



(17) 



+ 



_Ql 

(u, + §) 



2oj 
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Mum) 



(Ql/Ql 



In [|] a simple powerlike input distribution in k 2 (truncated at the scale of 
the limits of the evolution integration over k 2 ) is used, i.e. the inputs are 



ff\x,k 2 ) = Ax- e {k 2 ) {uJ0+1/2) Q(k 2 -Ql) (20) 

3?>(N,u,) = , (21) 

(N + e)(uJo - to) 

ff\x,k 2 ) = Ax- e {k 2 ) {uJ0+l/2) Q(k 2 -Ql)(l-Q(k 2 -Ql)) (22) 

?f\NM = A{Q ; \ (23) 

(N + e)(uJo - lo) 

The dimensionless input distribution has normalisation A (which carries di- 
mensions [— (co>o + 1/2)]), ujq < 1/2 ensuring that its momentum distribution dies 
off at large k 2 and e ~ 0.08 gives an x dependence motivated by the observed slow 



rise of the proton-antiproton total cross section at high energies (see e.g. [jig) ). 



With these choices of input distribution eqs. fll5jL7D agree with those derived in 



upto a minus sign (for the latter) as described in 
For later convenience we rewrite eq. (|T8"|) as 



N + e\ (w-a) i^ + V J 



where are Xx,X 2 are given by the terms in square brackets in eq.fllBD with the 
A/(N + e) factor removed. 

We now investigate the sensitivity of this model to the degeneracy of the upper 
and lower scales used to parameterise the input distribution and those employed 
to restrict the k 2 integrals in the evolution equation. Explicitly we replace Qq, Q\ 
in eqs. (|21|j23| ) with /Xo,/xf. By keeping these scales distinct and using the same 



input in each case we demonstrate the smooth matching between the solutions as 
Qq — > 0, Q\ —>■ oo. A comprehensive series of analytic and numerical cross checks 
are performed which give us confidence in the solutions. 
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However there is a more important reason for making this distinction. Physi- 
cally we expect the input distribution to be peaked around small values of k 2 with 
a small tail in the perturbative high-fc 2 region. It is these perturbative gluons 
which we may justifiable evolved in perturbative QCD using the BFKL equation. 
For convenience we choose to use a theta function to model the input, thereby 
restricting ourselves to the high-fc 2 gluons above some scale. Hitherto we have 
not made the distinction between this scale and the infra-red cutoff scale on the 
integral of the BFKL equation (which defines how the input distribution evolves 
with x). This is clearly an unphysical restriction since the choice of input should 
be independent of how we choose to model the kernel. 

We proceed now with the simple model in which these two scales are kept 
distinct. For generality, we now choose to use input J 7 ^' in the non-cutoff case 
too. We assume for the time being that /x 2 , > Qq since we are interested in taking 
the limit Ql -> 0. 

The o;-plane contour of eq.([3]) is closed in either the left or right half plane 
according to the behaviour demanded by the "closure factor" in the integrand. 
In the non-cutoff case the closure factor, (k 2 /^)^, forces closure to the left for 
k 2 > [Aq and gives 



A ( p(l/2+w ) p(l/2-u>jvi) 2(u> +awi)/l _ 2 ^ 



(N + e) \ (l + N^KiuJo)) 2lu n1 (lu + lu n1 ) 



(25) 



and to the right for k < llq 



A /rI-,,,2 \p(l/2+w m ), . 2 (">o-wjvi)\ 

UNk 2 ) =jj^— U T\ ^ • (26) 

{N + e) y 2u N1 {u - lo nx ) j 

In the single-cutoff case the factor (k 2 / nffi enforces closure to the left for the 
region k 2 > /xjj > Ql to give 



~ A ( k 2 ^) p(l/2-cJjvi)(j _ 

M ' J (N + e){ (l + N^K(u )) + 2u m 



I 2 + ^NlUh (.2 ~ u Nl)V0 (Voj 

(wo + ^ivi) (wq-^jvi) 



(27) 
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and to the right for /x 2 , > k 2 > Ql 

f s N,}t (N + e) 2uj ni {loq - um) 

((1 + uj m )k 2 ^ 2 ^ - (1 - w^)^ 172 -"* 1 ^^)*"* 1 ) • (28) 
Similarly the double cutoff case solution is given by 



UN,k 2 



A 



j^2(l/2+w ) 



+ 



X 



+ 



X 



(iV + e) \ (l + JV-^wo)) 

fc 2(l/2-„„ l)( l 



(UJ + U N1 ) 
^2(1/2+^0 (I 



A(a;jvi) 



2c; 



ATI 



(a> - ^ivi) 



A(aj m ) 



(29) 



for k 2 > n 2 > Q 2 , 



fd(N, k 2 ) 



(AT + e) 2^ 
+ 



(u; - ^JVl) 



/c 2(j r +^ivi)g-2a;iv 1 X 2 (uJ N i) fc2 (|-a; J v 1 )g2a; J v 1 ^"l(^JVl) \ 

1 A(cj Ari ) ' ~' J A(uj m ) / 



(30) 



for fil > k 2 > Ql, and 



fc 2 ) 



.4 



(iV + e) 
2uni 



^2(1/2-^0(1 _ w 2. J 



, 2(o)o+Wivi) 2(w +^jvi) 
/fO ~ A*l 

(^o + Um) 



+ 



■X 



f ^2(l/2-o;yi)g2a>jvi Xi(uJ N1 ) _ ^,2(l/2+o; wl )g-2a; wl X 2 (u! N l) \ 

\ A(a;jvi) ' A(a; JV i) / 



(31) 



for Qf > k 2 > {i\ > $ > Ql and X U X 2 are defined in eq.(|24|). f d (N,k 2 ) is 
the most general solution in that it uses the most general input distribution and 
contains the solutions in the other cases. 
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At first the solutions derived above seem very complicated. However one 
may check that they are indeed correct by looking at the solutions in various 
limits. With the definition fll2|) of uni as the positive square root we have that 
$le(u}Ni) > along the contour. We will use this result in what follows. First we 
look at / s in the distinct cutoff case. If we take /Xq = Qo i n eq.(|27|) we obtain 



Ak 2(l/2+U0) ^2(1/2^*0(1 

H ttt— — — — (32) 



;i + N- l K(u )) V (N + e) ' (N + e)(u + |) 

which is precisely what one would obtain upon inverting eq.(19) of ||]. In the 
limit Ql -> eqs.(|^|D tend to eqs.(H§|). So we have 



lim f,(N,k*,t4,Qt) 



hm f.iN^^Ql) 
Qi-*0 



f s (N,e,Qi) 



(33) 
(34) 



The double cutoff case is slightly more work but one may also prove 



\im f d (N,k 2 ,fil,Ql,f,l,Ql) — f d (N,k 2 ,Q 2 Q ,Ql) (35) 



hm UN,k 2 ^lQl^lQ\) _^ f a (N,k?,f4,Ql). (36) 

Q?->-ao 



Having seen the above cross checks working analytically we now have confi- 
dence in our solutions and proceed to invert the remaining iV-plane integral. To 
complete the inversion we need to perform the iV-space integral : 

f(x,k 2 ) = f ^-x N f(N,k 2 ). (37) 
Jc N 2m 

The contour, Cn, lies parallel to the imaginary axis and to the left of all sin- 
gularities in the iV-plane the nearest of which is the square root branch cut at 
iV = —K . Writing N = N R + zvVj, (with N R fixed along the contour, C N ) the 
x N factor may be written 

X N = x Nr [cos(Ni ln(x)) + % sin(vV/ ln(x))] 
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which tells us that the integrand is oscillating along the contour with period 
2tt/ lax. 

Let us now consider the symmetry properties of the integrand under the trans- 
formation Nj — ► — Nj. For any analytic function / which is real on some part of 
the real axis (which is true for the functions we will be considering) we have by 
the Schwartz reflection principle: 

f(N R + tNj) = MeifiNR + iN^+iQmifiNR + iNj)) 
^e(f(N R -iN z )) = Ue(f(N R + iNj)) 
ZmifiNR-iNj)) = -^m(f(N R + tNj)). 

We see that the imaginary part of the integrand in eq. (|37|) is odd overall and 
the contributions from above and below the real axis cancel to produce a purely 
real function for the gluon density. So we have 



f(x,k 2 ) = 2x NR — - cos(AT 7 lnx)Ke(/) - sin(JVjlnx)9m(/) (38) 

JO 2lX l J 

The non-trivial iV-plane singularity structure associated with the square root 
branch points at N = — K and N = contained in the factor 1/2o;jvi necessitates 
a numerical inversion of the remaining transform. We achieve this using a Fortran 
program which calls the NAG integration routine DO 1 ASF which is specifically 
designed for integrals with oscillating integrands. Having performed the iV-plane 
integral we have explicitly checked the smooth matching conditions for all ranges 
of k 2 and x. 

The small-x behaviour, in the non-cutoff and infra-red cases, is dominated by 
the square root branch point at iV = — K (c.f. eq.(|T2|) ). In the double cutoff 
case the singularity structure in the iV-plane is governed by A(um) = defined 
in eq . (|19|) , the zeros of this function occur at pure imaginary values of ujni and 
also at ujni = [|], [Tj. In H it was shown that for the double-cutoff case this 
branch point singularity cancels between terms to leave a softer behaviour given 
by a sum of poles to the right of K in the iV-plane. 

It appears from eq.(|29|) that we now have a singularity in fd{N, k 2 ) at iV = 
— K associated with the l/2um factor. A careful consideration of the terms in 
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square brackets, in the limit uni — ► 0, however reveals that this singularity is 
cancelled and one is left with the sum of poles in the iV-plane coming from the 
zeros of A(o>jvi) for which um are imaginary as before ||, [/J. 



3 Using the full BFKL kernel 

In this section we explain how to implement a more realistic model which uses the 
analytic form for the full BFKL kernel given by eq.([8|) (K(u) = a s x(uj)). Initially 
we include only the nearest poles at ±1/2 in the kernel. These are referred to 
as the leading twist poles since, as we will see, they lead to the largest power in 
Q 2 . The pole positions in the tu-plane, at ±o?jvi, are given by the solution to the 
following simultaneous equations 



N R + Ue(K(u m )) = 

Nj + $sm(K(u> Nl )) = (39) 

such that 

\3te("Ni)\<\ (40) 

The exact position of these poles depends on the choice of position of the contour 
in the iV-plane, and on the position along this contour. Looking at the graph 
of $le(K(cj)) in fig.(|l|) we see that the solutions of the real part of the equation 
are given by a curves of constant height, —N R , such that — Nr > K Q , where 
Kq is the height of the saddle point at u — (0, 0) (imagine taking a horizontal 
slice through the saddle). Provided this condition is satisfied, N R lies to the 
left of all singularities in the iV-plane. As Nj varies the solutions ±u>ni move 
along these curves. Having chosen the position of the iV-plane contour we find 
cjjv! numerically, the complex solution then feeds into the residue of the poles 
calculated below. 

We may now invert the solutions ([7,15, 17) to (N, fc 2 )-space, by closing the 



contour in the appropriate direction and summing over the residues of the enclosed 
poles as before. To do this we need to determine the contribution that the crucial 
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factor, 1 + iV _1 i^(ci;), makes to the residue. For example, if a particular term 
contains a closure factor forcing closure to the left we need to know the residue 
of the pole at u — —ujjsiii so we expand the factor about this pole to arrive at 



1 



1Z(u) 



1 + N- x K(u) 
U(-u m ) 



(<jj + L0 NX ) 

d 



£( 1+ >>) 



(a s /N) x '(-u>m) 
duj 



UJ = — LJJVl 



(41) 



with similar relations for the pole at u — ujni- The function xi 00 ) which retains 
the full cu-dependence of the kernel is given by eq.@ 

For the non-cutoff case, closing to the left for k 2 > Hq, we have 



UN, k 2 



AN 



£2(l/2+wo) j^.2(l/2-w m )^2(wo+wjVl) 



(N + e) 
and to the right for k 2 < 

fn(N,k 2 ) 



+ 



(N + K(uj )) (uo + um)oc s x'(-um) 



(42) 



^2(1/2+^1)^2(^0-^1) 



(N + e) y [Uo - . • • \ 
With an infra-red cutoff on the evolution, Qq, we get 



(43) 



L(N,k 2 



AN 
(N + e) 



£,2(l/2+cj ) 



j^.2(l/2-w m ) 
+ rX 



(AT + K(w )) a sX '(-u m ) 

2(w +^jvi) / 1 ., \ , ,2(o)o— Wjvi)/^2(2wjvi) 
Mo 1 2 ~~ u Nl)^0 vo 



(44) 



(cJq+CJjvi) (| + ^jvi)(^o-^jvi) 

for k 2 > fiQ > Qq. For /Iq > k 2 > Ql the input term is closed to the right to give 

fs(N } k 2 ) 



AN f^-^) 



x 



(N + e) a sX '(-um) 

fc2d/2+^i) _ (j_ a;jvl ) fc 2(l/2-a, Jvl )g2(2^ 1 ) ' 

(a»o-^jvi) (| + ^tvi)(^o - ^jvi) 

Similarly the double cutoff case solution is given by 

AN ( fc2(l/2+c^o) J 



(45) 



+ 



(N + e)\ (iV + e)(iV + K(u; )) a.^-u^) 
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X 



+ 



2(wo+wjvi)^2(l/2-w is ri) ^2(a;o-w A ri)^2(l/2+£J i vi) 



(u + u N1 ) (uq-vni) 
a s x'(-uJ N1 )^Ni)(h + uni) 



(46) 



for Q\ > fil > k 2 > $ > Ql, 

AN 



fd(N, V) 



+ 



(N + e)a s x'(-uJ N1 ) y (uj q - u N1 ) 



(47) 



for Ql > fil > $ > k 2 > Ql, and 
UN,k 2 ) AN 



'^.2(l/2-wjvi)^2(wo+Wivi 



2(tJo+wjvi)\ 



(iV + e)a s x'(-^ivi) 



(U + Wjvi) 



+ 



A(u;Ari)(u;Ari + |) 



(4* 



for Q 2 > k 2 > fi 2 > fi 2 > Q 2 . 

Note that in these equations only terms corresponding to the input distribu- 
tion change in the different fc 2 -regimes. So the first and second terms in eqs.(H3, 
46]) are the same as eq.( fl^) ; the first terms in eqs. ([45| , are the same as eq.([43|). 
The third term in eq.([46]) and the second in eq . fl4~7|) (i.e. those involving /j, 2 ) are 



present because a different input is used in the double cutoff case. The remaining 
terms in eqs.(i"l, |4l)| , |4"6] , |4"T| ) represent the effect of the cutoffs on the evolution. 

The remaining iV-plane integral is now performed numerically as described 
above employing extra subroutines in the Fortran programs to calculate K{uq) 
and x\ u ni)- Again we verify that the numerical solutions have the expected 
properties, i.e. that they satisfy the various limits detailed above. This has been 
done for k 2 > /Xq and k 2 < /Jq corresponding to the inversion of eqs.(42, 44, 46) 
and eqs.(|43|, |45| , |47D respectively and we find that all the appropriate cross checks 
are observed. 
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Fig. (g) shows the effective small- x slope, A e g, defined by 

_ d In G(x,Q 2 ) 
XcS = d\nl/x ■ (49) 

in each percentage of its asymptotic value of Ao = 41n2a s . Two curves 

are shown for the double cutoff case corresponding to the range of variables 
expected from HERA kinematics. The two values chosen for Q\, 10 3 and 10 5 
GeV 2 , correspond (somewhat arbitrarily) to the maximum value of Q 2 at x = 10~ 2 
and the total centre-of-mass energy respectively. 

3.1 Higher twist terms in the leading ln(l/x) approxima- 
tion. 

We now consider including further poles associated with the BFKL kernel. As 
explained above these poles are subleading in Q 2 and may be referred to as "higher 
twist" terms. Since the solutions get very complicated with more poles included, 
we restrict ourselves to the non-cutoff and infra-red cutoff cases as these contain 
the essential features. 

Keeping the nearest and next-to-nearest poles in the kernel (w = ±1/2,3/2) 
leads to four solutions of eq.(|H]) at uj = ±Uni, ±wjv2 where 

1 < \Ke(u N2 )\ < | (50) 

this can easily seen from fig.(|l]) which shows the full kernel out to uj = ±3/2. We 
note that 3?e(u;jv2) ~ 3?e(o;jvi) + 1. 

The essential features can be illustrated using the following simple four pole 
kernel which has the same pole structure as the full kernel close to the contour: 

^ = (T^) + (9^)' < 51 > 

where K\, K 2 are determined, for example, by matching this kernel onto the full 
kernel (eq.(||)) at uj = 0, ±1. Eq.([T0|) now has four solutions given by the solutions 
to a quartic equation in ujn- The residue factor at these poles is given by (c.f. 
eq.@) 
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Figure 2: The effective slope of G(x, Q 2 ) as a percentage of its asymptotic 
value A = 41n2a s , for Q 2 = 10 GeV 2 . The solid line in each pair of curves 
is for Ld = —1.0 and the dashed line for uj = —2.0. At small enough x the 
slopes approach a value independent of the specific choice of cu . For the curves 
shown in the single cutoff case (labelled V),Qq = 0-5 GeV 2 . In the double 
cutoff case we show two sets of curves; for d\, Q\ = 10 5 GeV 2 whereas for d2, 
Q\ = 10 3 GeV 2 . fx 2 is set very high (10 9 GeV 2 ) so that the comparison can be 
made for the same input distribution (here /Xq = 0.5 GeV 2 ). 
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U(uj n ) 



l)(4^-9) 



(52) 



8lu n [K^l 



9) + K 2 {4u 2 N -l)\ 



which becomes singular either when uo^ = or when the term in square brackets 
in the denominator becomes zero. For the leading poles at uo = ±uo^i the former 
condition is satisfied when N = —Kq = — 41n2a s i.e at the original branch point. 
For the subleading poles in the uo -plane at uon = ±^tv2 the latter condition may 
be satisfied, this occurs for two complex conjugate values of N rs — iT (0.1±0.25). 
Since these singularities lie far to the right of the branch point in the iV-plane at 
iV = —Kq one expects them to be insignificant for small enough x. 

This analysis may be extended trivially to include further poles in the kernel. 
The pole structure of the kernel determines that the solutions to iV + K(uni) = 
are spaced approximately one unit apart in terms of dte(uNi). This leads, on 
closing the uo plane contour, to a sum of k 2 terms with each subsequent term being 
approximately one power in k 2 down on the previous one (or equivalently in Q 2 
once the k 2 integral has been performed). In this sense all the terms after the first 
one are higher twist in the language of the operator product expansion and may 
be neglected for large enough values of Q 2 . In what follows we investigate, and 
attempt to quantify, the magnitude of the contribution made by these additional 
poles at intermediate values of Q 2 and x, by explicitly including four poles in the 
BFKL kernel. This will be referred to as the next-to-leading pole calculation and 
denoted with a superscript "nl" (with a superscript "1" now for the leading pole 
solutions). 

The non-cutoff case solution is straightforward. The effect of including this 
pole is to introduce terms in eqs. (|4~2"1J4"3D in uo^2 identical in form to the existing 
terms in u^i, which are, approximately speaking, of order k 2 / '/Xg relative to the 
leading, k 2 ^ 1/2 ~ UJN1 \ terms. 

Following the procedure outlined in [/]], in the infra-red cutoff case the function 
S S (N, uo) is determined by the fact that is removes the right half plane poles (at 
uo = u}ni,un2) from T s . Assuming a form which has the appropriate closure 
factor 
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one solves the following simultaneous equations: 



J*°\N,u N1 )+S(N,u m ) = 
^°\N,u N2 )+S(N,u m ) = 



to give 

(1/2-^0(1/2-^2) 

x '(3/2 -^og^ 1 ^^^) 

(3/2-^0(3/2-^) 

(io N2 — ujm) 
Xl/2-u m )Ql^T®{N,u m ) 



x 



{3/2-u N2 )Ql^^ \N,u N2 ) 



[1/2-lu N2 )QI"™FW(N,u N2 ) 



(53) 



We can now see that both a>\ and a 2 contain information about both the leading 
and next-to-leading poles. As we will see, this has a significant effect on the x 
dependence of the solution. 

After inverting the cu-plane transform as appropriate for the values of k 2 (using 
the full kernel, K(uS) = a 8 x(u), of eq.(§)) we arrive at the following solutions: 

AN ( fc 2(1/2+u,o) 



f s (N,k 2 



(N + e) V 

, 2(wo+wjvi) 
P0 



£,2(l/2-u) m ) 

(N + K(u )) + a sX '(-um) X 



(C0 + Wjvi) 



+ 



«2 



(1/2 + ^x) (3/2+Wjvi) 



(54) 



for k 2 > fil > Qq. For fi^ > k 2 > Ql we have 



MN, k 2 



AN 



x 



^.2(l/2+wjvi) ^2(aJo-wjvi) 



^2(1/2-0)^) 



-g 



2o;jvi 



a 1 



+ 



«2 



(1/2 + ^0 (3/2+wjvi) 



(55) 



Until now we have been working with the unintegrated gluon distribution 
function f(x, k 2 ). This is not a physical observable. It must be convoluted with 



an appropriate coefficient function to arrive at an observable for the particular 
process under consideration. For simplicity, we present our results in terms of 
the gluon structure function G(x, Q 2 ) defined in eq.([|). Since we have analytic 
expressions for f(N,k 2 ) in each case we may perform this k 2 integral prior to 
inverting the iV-plane integral. 

In the non-cutoff case we split the integral up into two regions and integrate 
eqs.(H,|3D, 

G n (N, Q 2 ) = f ^-f n (N, k 2 ) + f ° ^-f n (N, k 2 ) (56) 



4 k 2 jnK ' 1 Jo k 2 



which gives for the leading pole solution 

AN ( (Q 2 (V2+-o) _ ffiW+wh 



G l n (N, Q 2 



+ 



(N + e) y (lu + 1/2)(N + K(lu )) 

2(1/2+cj ) / -, 



asX'(-um) \(l/2 + wjvi)(w - am) (1/2 - uj n1 )(uj + u N1 ) ^ 



d s X'(-uJNi) I (1/2 -u N1 )(vo +um) 
The next-to-leading solution has an additional term relative to this 



AN 

G n n l (N,Q 2 ) = G l n (N,Q 2 ) + -^—^x 



2(l/2+u ) / i -, 



oisX'i-um) \(l/2 + u m )(u - oj N2 ) (l/2-Lo m ){uj Q J ru m ) j 

Ot s X'{-^N2) y (1/2 - UJ N2 ){u Q + UJ N2 ) ) J ' 

In the single-cutoff case, which has the same input, the integral is non-zero in 
two regions: 

Gs(N, Q 2 ) = /J ^f n (N, k 2 ) + J* ^f n (N, k 2 ) (59) 

the result differs from the non-cutoff solution only in the coefficient of Q 2 ( 1 / 2 -^ni) 

AN I Q 2 (V2-ivi) M 2 ^-^)(g 2 ) 2 ^ 



G l s (N,Q 2 ) = G l n (N,Q 2 ) - 



(N + e) ya s x'(-u N i) (1/2 + u m )(u - w m ) 

(60) 
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The next-to-leading solution can no longer be expressed in terms of the leading 
solution plus an extra piece from the new pole as in eq.(^8[) as a result of the 
information about the position of both poles tied up in the coefficients a\, a 2 . The 
solution is given by 



G n s l (N,Q< 



G n n l (N,Q 2 



(Q 



2(1/2-^1) _ q2(1/2-w N i) 



(1/2 -u N1 )a s x'(-wm) 



1/2 + um 3/2 + ^1 



(61) 



For completeness we also give the double cutoff solution for the leading pole 
case for the same fc 2 -regimes as in the single cutoff case: 



G l d (N,Q 2 ) = G l n (N,Q 2 ) + 



AN 



Q 



2(l/2+w m) 



. 2(o>o— o)jvi) 



+ 



/ - J 2(l/2+w^ri) / 2(ojo— wjvi) ,, 2 (wo— 

Vo (Po ~ n )_ 

a s X'(-u} N1 )(l/2 + uj n1 )(ujo - u N1 ) 
1 

a s X / (-o; i vi)(l/2 + u N1 )A(u N1 ) 

( g2(i/2-^ 1 ) _ gg(i / 2 -^))g^ix 1 (cu jV1 ) 
(1/2-wjvi) 
{Q 2(l/2+» N1 ) _ Qf/^—))Q~^X 2 {u m ) 



(62) 



(1/2 + ^) 

We now compare the model with two poles in the kernel to that with four in 
order to determine the significance of these next-to-leading poles as a function 
of Q 2 and x. To this end we work with the following quantity (which is the 
percentage difference between the two solutions) 



100 x 



GW(x,Q 2 )-G«\x,Q 2 ) 



G«(x,Q 2 ) 

From the above analysis we expect R n to go to zero as Q 2 
coefficients of the leading Q 2 piece (the term in g 2 ( 1 / 2 -^Jvi)) are the same in each 



(63) 

00 since the 
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case (see eqs.( p)T||)%D ). In contrast, the coefficients of the "leading twist" pieces 
for the single cutoff case are different (see eqs.(^,^l|)) so we expect R s to tend to 
a constant (which will be determined by the choice of input parameters u>o, //q, Qq] 
the overall normalization cancels in the ratio). Figs.(||,[|) illustrate these features 
for a specific choice of input parameters. 

Let us now consider the x dependence of these solutions. We have seen from 
the simple four pole model of the kernel that the singularities in the iV-plane asso- 
ciated with un2 are sub leading compared to the branch point at iV = — Kq. The 
non-cutoff solution contains terms which involve either uj^i or but not both, 
hence one expects that the significance of the next-to-leading poles will decrease 
as x — > 0. One can see from fig.(|5|) that this is indeed the case. In the single 
cutoff case however the effects of the poles are coupled in the coefficients a±,a2, 
this makes the x dependence of R s considerably less trivial since the coefficient 
of the leading piece in x in and G® are different. 

Fig. ([]) shows R s as a function of x, for various Q 2 and two different values of 
the input parameter ujq. uq = — 1.0 is special in that it lies in the same interval 
in u as co>tv2- The curves shown for u = —2.0 are typical of those for which u 
lies outside of this interval (< —1.5). For the latter, R s is seen to increase as x 
decreases eventually tending to a constant percentage (of the order of 10 %) as 
x -> 0. 

It seems that for intermediate values of x when an infra-red cutoff is present 
the next-to-leading poles in the BFKL equation are significant and including only 
one of the subleading poles leads to a behaviour in x similar to what one would 
expect for the sum of two powers. 

This conclusion follows given the fact that we have chosen to regulate the 
infra-red physics using a theta function cutoff on the integral in eq.(0) and a 
specific form fl20D for the input distribution. It would certainly be interesting to 
establish whether it is true in general. 
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Figure 3: R n as a function of Q 2 for various values of x showing that R n — > 
as Q 2 — > oo. Here = 0.5 GeV 2 and cj = —2.0 . 

Conclusions 

We have presented analytic solutions to the BFKL equation with infra-red and 
ultra-violet cutoffs in moment space, including the leading pole of the full BFKL 
kernel. We have investigated the numerical significance of the sub leading ('higher- 
twist') poles. We find that when one includes an infra-red cutoff the assumption 
that these poles are irrelevant is brought into question. 
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Figure 4: R s as a function of Q 2 for various values of x. R s tends to a constant 
as Q 2 -> cx). Here ^ = q.5 GeV 2 , Q 2 , = .5 GeV 2 and cu = -2.0. 
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Figure 5: R n as a function of x for various values of Q 2 . The percentage 
difference between the two solutions is less than 2% and tends to zero as 
x — > 0. The curves shown correspond to (largest first) Q 2 = 10, 100, 1000 
GeV 2 . Here $ = 0.5 GeV 2 . 
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Figure 6: R s as a function of x for various values of Q 2 . For the solid curves 
lu = —2.0, wheras for the dashed curves u = —1.0. In each case the curves 
correspond to (largest first) Q 2 = 1, 10, 100 GeV 2 . The input parameters are 
ii\ = 0.5 GeV 2 ,Q 2 = 0.5 GeV 2 . 
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